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ABSTRACT 

Wc present a phenomcnological study of highly ionised, non-photospheric absorption 
features in high spectral resolution vacuum ultraviolet spectra of 23 hot DA white 
dwarfs. Prior to this study, four of the survey objects (Feige 24, REJ 0457-281, G191- 
B2B and REJ 1614-085) were known to possess these features. We find four new 
objects with multiple components in one or more of the principal resonance lines: 
REJ 1738+665, Ton 021, REJ 0558-373 and WD 2218+706. A fifth object, REJ 2156- 
546 also shows some evidence of multiple components, though further observations 
are required to confirm the detection. We discuss possible origins for these features 
including ionisation of the local interstellar environment, the presence of material 
inside the gravitational well of the white dwarf, mass loss in a stellar wind, and the 
existence of material in an ancient planetary nebula around the star. We propose 
ionisation of the local interstellar medium as the origin of these features in G191-B2B 
and REJ 1738+665, and demonstrate the need for higher resolution spectroscopy of 
the sample, to detect multiple ISM velocity components and to identify circumstellar 
features which may lie close to the photospheric velocity. 

Key words: stars: abundances - stars:atmosphcrcs - stars:white dwarfs - circumstel- 
lar matter - ultraviolet:stars 



1 INTRODUCTION 

Spectral lines from highly ionised species are observed 
at non-photospheric velocities in several hot DO white 
d warfs. Of the 11 star s included in the survey presented 
bv lHolberg et all il998T) (hereafter referred to as HBS), six 
objects, all with temperatures T e fr ^ 70,000 K, were found 
to ex hibit highly ion i sed fe atures at non-photospheric veloc- 
ities. iHolberg et alJ il998ft argue that such features cannot 
be of an interstellar origin since such highly ionised species 
are uncharacteristic of the local ISM, and are not observed 
along adjacent lines of sight to stars at greater distances. 
Further, while interstellar lines can be observed at velocities 
which are blue- or red-shifted with respect to the photo- 
spheric value Wphot, only blue-shifted features are found in 
the DO sample discussed by HBS. Such features provide 
strong evidence for the existence of ongoing mass loss in hot 
DO white dwarfs. 

Although over 80% of know n white dwarfs belong to 
the DA group JSion et al.lll997fi . relatively few have been 
found to exhibit non-photospheric features. Only 5 of the 44 
DA stars considered in HBS are accompanied by any type 
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of circumstellar feature, with no apparent dependence on 
temperature. In some instances, the appearance of this phe- 
nomenon can be plausibly explained in terms of interactions 
between binary components, as in the DA + dM 1.5 sys- 
tem Feige 24, which shows multiple components in the CIV 
AA 1548, 1550 doublet. 

However, in the DA white dwarf CD- 38° 10980, which 
is a member of a wide binary system, IHolberg' et al.l (I1995D 
find Si and C absorption features in IUE spectra, which 
are shifted by -12 km s _1 with respect to the photospheric 
velocity of the with dwarf, which can be inferred from its 
measured gravitational redshift. Holberg et al. showed that 
the atmosphere of this object was devoid of Si and C, at 
the expected photospheric velocity, and used the presence of 
excited or meta-stable levels in the shifted features as evi- 
dence that the material was not located in the ISM along the 
line-of-sight to the star. These observations were explained 
in terms of a dense, gaseous halo in close proximity to the 
star, possibly an extension to the atmosphere, for which a 
similar temperatur e and e lectron density was derived. Alter- 
nately, IWolff et alJ fcOOll) using FUSE observations of CD- 
38° 10980, observed a set of Si III lines shortward of 1120 A 
which they attribute to the stellar photosphere. They were 
able to successfully model the observed equivalent widths 
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of these lines as well as the Si lines seen in the IUE obser- 
vations with a photospheric Si abundance of 2xl0 -8 . How- 
ever, the velocity discrepancy remains unexplained. In an- 
other example of circumste llar features in an isolated DA 
object. Irlolberg et alJ il997Tl found evidence for weakly blue- 
shifted CIV and Si IV components in REJ 1614-085 (T eff ~ 
38,500 K) at ~ 30% of the strength of the photospheric lines, 
shifted by -25 and -40 kms -1 respectively. Similar features 
in the spectrum of the T efr ~ 57,000 K DA, RE J 0457-289 
have also been discussed bv lHolberg et all lll997l) . 

Agreement between the predicted and observed abun- 
dances of atomic species in the atmospheres of white 
dwarf stars has improved with the introduction of strat- 
ified model atmosphere codes, beginni ng with the strati- 
ficatio n of He and Fe investigat e d bv iBarstow fe Hubenvl 
]l99Sf) and IBarstow et alJ il99Sft . iDreizler fe Schuhl J200lf> 
describe models in which stratified abundances are calcu- 
lated self-consistently, by considering the depth dependence 
of temperature, density, radiation field and level populations 
using an iterative procedure. These models have succeeded 
in reproducing the soft X-ray, EUV and FUV spectra of 
a sample of DA stars, and can explain the widely varying 
levels of metallicity in hot DA white dwarfs. However, this 
agreement is not complete; for example, lines of C, N and 
O are not reproduced as accurately as those of Fe and Ni, 
while observed differences b etween objects of simila r T e g 
and log g are unexplained iSchuh fc Dreizleill200lf) . The 
success of stratified models is encouraging, but the obser- 
vation of highly ionised circumstellar features suggests the 
existence of processes which may modify the predicted equi- 
librium abundances, and improving the agreement between 
observed and predicted abundances requires better under- 
standing of the nature of these features. It is particularly 
interesting that the DO white dwarfs, in which circumstel- 
lar features are relatively common, are al so poorly modeled 
by th e new generation of stratified codes (IDreizler fc Schuhl 

Data from the STIS instrument on board HST have al- 
lowed the study of white dwarf spectra to be carried out with 
an accuracy impossible to achieve using earlier instruments, 
and at a resolution which permits more precise examination 
of intrinsic line profiles. A more sensitive search can now be 
made for signs of mass loss and accretion which may modify 
equilibrium abundances in white dwarf envelopes. 

In section we present a phenomenological study of a 
sample of 23 hot (T cff ^ 20,000 K) DA white dwarf stars 
for which either GHRS, STIS or high resolution IUE echelle 
spectra were available. In each star, the resonance doublets 
of CIV, NV and Si IV have been examined for signs of mul- 
tiplicity (whether in the form of asymmetry or distinct com- 
ponents), and statistical tests applied to determine the sig- 
nificance of proposed secondary features. In several stars, 
Gaussian profiles are used to model the observed lines. A 
Gaussian approximation to line data can be justified since, 
at the resolution of both the IUE echelle data (R~20,000) 
and the STIS E140M grating (R~40,000), unsaturated ISM 
lines are not resolved, nor are the photospheric profiles of 
most absorption lines. In some cases, elemental abundances 
are considered as additional evidence for the presence of non- 
equilibrium processes. The discussion of individual objects 
is divided into two subsections, covering stars which show 
clear (or suspected) highly ionised, non-photospheric fea- 



tures, followed by those that do not at the resolution of the 
data used in this study. Several of the stars in this category 
exhibit unusual features which require further investigation 
at higher spectral resolution. 

In section the results of this study are discussed, and 
a variety of explanations for the presence of highly ionised 
non-photospheric features in certain white dwarfs, are con- 
sidered. 



2 OBSERVATIONS AND ANALYSIS 

2.1 Observational data 

T he sample o f star s was chosen to match that 
of lBarstow et al] i200ll) . with the addition of the super-hot 
DA PG 0948+534 and Ton 021, since GHRS, STIS or 
high resolution IUE echelle spectra were available for 
these objects. Table summarises the survey stars, their 
basic physical parameters, and the source(s) of data used 
in this wor k. Values for t e mpera ture and gravity were 
taken from IBarstow et al] l [200lt). The adopte d visual 
magnitudes were those of iMarsh et al] il997f) . except 
where otherwise stated. The mass, radius and distance 
to each star was estimated using the evolutionary models 
developed by IWoodl dl995h . taking the stated values of 
T c g , log g and m v as input parameters. 

The results of the survey are summarised in ta- 
ble [5] which includes the measured velocities of interstel- 
lar («ism), photospheric (w p hot), and any non-photospheric 
highly ionised lines {v cil - c ) for all stars in the sample. Note 
that values for v circ are not relative to the photospheric fea- 
tures, but are absolute velocities. Also included in the table 
are estimated values for the escape velocity (v BBC ) and grav- 
itational redshift (v srav ) of each star, and the velocity of the 
primary component of the local interstellar cloud along the 
line of sight to each star (f L ic). As noted previously, the 
spectral resolution of IUE, GHRS and STIS (in the E140M 
configuration prevalent in this study) is insufficient to com- 
pletely resolve the ISM components, and hence the value of 
vism presented in table [5] represents the velocity of the pri- 
mary component (or blend) observed in the data. This table 
also includes calculated values for gravitational redshift, es- 
cape velocity, and the velocity of the principal component 
of the local interstellar cloud in that direction. The local in- 
terstellar cloud can be described with reasonable accuracy 
by a cloud moving at 26 ± 1 kms -1 (heliocentric velocity) 
towards /„ = (186 ± 3)°, bu = ( -16 ± 3)°, or a = 74.5°, 
S — +15° llLallement et alJll995f) . The velocity of absorp- 
tion lines from the LIC in a particular star, i> L ic*, may be 
estimated from the projection of u L ic onto the target direc- 
tion. 

2.2 Analysis methods 

Absorption line parameters were measured using an IDL 
code, "Lines", written by one of us (JBH). The routine mea- 
sures the properties (wavelength, velocity, equivalent width 
and associated uncertainties) of cursor-defined features in an 
input spectrum, given a user-supplied rest-frame wavelength 
for the feature. The principal lines in each spectrum were 
identified using the data contained in HBS, beginning with 
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Table 1. Stars includ ed in the surv ey of circumstellar features, in descending tempera ture order. Mass, radius and distance calculated from the 
evolutionary models of lWoodl il995f) . Visual magnitudes taken from lMarsh et all Jl997tl unless otherwise stated. 



Star 


en \ j 


1o£ 


M (Ma,) 


R (Rm) 


ttl v 


D (pc) 


Data source [Mode] 


Resolution 


PG 0948+534 


90,000 


7.27 


0.478 


0.027 


13.71 e 


193.8 


STIS [E140M] 


40000 


RE J 1738+665 


71,300 


7.53 


0.540 


0.024 


14.61 


243 


STIS [E140M] 


40000 


Ton 021 


69,700 


7.47 


0.550 


0.023 


14.539 


217 


STIS [B140M] 


40000 


REJ 0558-373 


63,000 


7.66 


0.580 


0.019 


14.37 


295 


STIS [B140M] 


40000 


REJ 2214-492 


62,100 


7.23 


0.464 


0.027 


11.71 


69 


IUE [swp] 


20000 


REJ 0623-371 


59,700 


7.00 


0.399 


0.033 


12.09 


97 


IUE [swp] 


20000 


WD 2218+706 


56,900 


7.00 


0.396 


0.033 


15.40 c 


436 


STIS [E140M] 


40000 


Feige 24 


56,400 


7.36 


0.546 


0.021 


12.56 


78 


STIS [B140M] 


40000 


REJ 2334-471 


54,600 


7.58 


0.536 


0.020 


13.44 


104 


IUE [swp] 


20000 


G191-B2B 


54,000 


7.39 


0.510 


0.021 


11.73 


50 


STIS [E140M] 


40000 


GD 246 


53,700 


7.74 


0.646 


0.016 


13.09 


72 


STIS [E140M]/IUE [SWP] 


40000/20000 


REJ 0457-281 


53,600 


7.80 


0.612 


0.016 


13.95 


108 


IUE [swp] 


20000 


PG 1123+189 


52,700 


7.52 


0.511 


0.021 


14.13 


117 


STIS [E140H] 


110000 


HZ 43 


49,000 


7.90 


0.561 


0.018 


12.99 


71 


IUE [swp] 


20000 


REJ 1032+532 


46,300 


7.78 


0.585 


0.016 


14.46 


127 


STIS [E140M] 


40000 


REJ 2156-546 


45,900 


7.74 


0.568 


0.017 


14.44 


129 


STIS [E140M] 


40000 


PG 1057+719 


39,555 


7.66 


0.519 


0.018 


14. 68 6 


411 


GHRS [G160M] 


22000 


REJ 1614-085 


38,500 


7.85 


0.596 


0.015 


14.01 


86 


GHRS [G160M] 


22000 


GD 394 


38,400 


7.84 


0.592 


0.015 


13.08 


57 


IUE [swp]/GHRS [G160M] 


20000/22000 


GD 153 


37,900 


7.70 


0.520 


0.017 


13.35 a 


73 


IUE [swp] 


20000 


GD 659 


35,300 


8.00 


0.662 


0.013 


13.37 


53 


STIS [E140M] /IUE [SWP] 


40000/20000 


EG 102 


20,200 


7.90 


0.573 


0.014 


12.76 d 


25 


IUE [swp] 


20000 


Wolf 1346 


20,000 


7.90 


0.572 


0.014 


11.52^ 


li 


IUE [swp] 


20000 



Bohlin et al J <1995h 
GreensteinUl984h 
Saffer et al.1 <1998ft 



Holber^^^al 



Holberg et al 



1997 ) 
1995ft 



c Naoiwotzki { 
f iDahn^t^lT 



jj ichonbernerl ^ 



Table 2. Measured velocities for interstellar, photospheric and non-photospheric lines for surveyed stars. Tentative identifications are 
marked "?". All velocities in kms - . 



Star 


V<,BC 


Ugrav 


Vlic 


VlSM 


«phot 


^circ 


PG 0948+534 


2645 


11.3 


10.00 


-0.26 ±1.26 


-14.25 ± 0.22 




REJ 1738+665 


3365 


14.2 


-3.46 


-18.56 ±0.09 


30.19 ±0.27 


-17.80 ±0.33 


Ton 021 


3074 


15.2 


10.61 


0.86 ±0.01 


37.14 ± 0.21 


7.62 ±1.27 


REJ 0558-373 


3477 


19.6 


25.19 


11.61 ±1.40 


22.74 ±2.81 


7.15 ± 0.77 


REJ 2214-492 


2526 


10.9 


-8.54 


-1.72 ±0.51 


33.49 ±0.45 




REJ 0623-371 


2143 


7.7 


14.46 


16.40 ±0.70 


41.15 ±0.56 




WD 2218+706 


2143 


7.6 


4.93 


-10.11 ±0.19 


-38.69 ±0.23 


-16.34 ±0.67 


Feige 24 


2588 


17.0 


20.52 


8.22 ±0.10 


80.33 ±0.50 


AV=25.0 ±3.5 


REJ 2334-471 


3253 


17.0 


-2.62 


15.15 ± 1.52 


40.57 ± 2.07 




G191-B2B* 


2679 


15.4 


20.58 


8.6 


25.20 ±0.26 


7.56 ±0.19 


GD 246 


3498 


25.7 


2.40 


-5.78 ±0.12 


-13.29 ±0.25 




REJ 0457-281 


3749 


24.3 


18.92 


14.48 ±2.15 


76.91 ±0.83 


21.76 ±1.27 


PG 1123+189 


3111 


15.5 


-0.36 


-0.67 ±0.06 


12.55 ±0.53 




HZ 43 


4461 


19.8 


-8.83 


-16.86 ±2.07 






REJ 1032+532 


3663 


23.3 


7.41 


0.84 ±0.21 


38.16 ±0.40 




REJ 2156-546 


3606 


21.2 


-9.65 


-8.39 ±0.18 


-19.94 ±0.68 


-1.65 ±0.76 ? 


PG 1057+719 


3384 


18.3 


6.51 


-2.89 ±0.69 


75.35 ±2.59 




REJ 1614-085 


3845 


25.3 


-25.29 


-29.56 ±0.33 


-37.31 ±0.40 


-71.42 ±6.63 


GD 394 


3801 


25.1 


-2.16 


-7.28 ±1.42 


28.75 ±0.91 




GD 153 


3444 


19.4 


-8.65 


-8.42 ±2.68 


12.45 ±2.07 ? 




GD 659 


4254 


32.4 


6.12 


9.77 ±0.22 


34.28 ±0.17 


-2.97 ±3.00 ? 


EG 102 


3934 


26.0 


1.05 


-1.37 ±3.85 


-0.69 ±1.44 




Wolf 1346 


3934 


26.0 


-10.76 


-15.38 ±0.95 


24.32 ±1.41 





*ISM velocity from ISahu et al.1 Il999t) : no error quoted. Value for ISM velocity determined from 
current work (using lower resolution data) is Uism = 16 ± 1 kms -1 . 
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unambiguous features such as the saturated interstellar lines 
of, e.g., Nl, and the resonance doublets of photospheric CIV 
and Si IV. Subsequent identifications were validated against 
the resulting ISM or photospheric velocity. 

Gaussian line profiles have been fitted to observed fea- 
tures to determine the velocity of circumstellar material, 
or multiple ISM clouds along the line-of-sight. The profiles 
were fitted using further bespoke IDL routines which applied 
first a single, and then dual, Gaussian to cursor- identified 
features in the spectra via the \ 2 minimisation technique, 
generating values for the velocity and equivalent width of 
the best-fit components. In cases where the dual-Gaussian 
fit was not obviously superior, x 2 values from the single 
and dual fits were then used to p erform a stand a rd F- test, 
following the method outlined bv lHolbere et all l|l 9971) . to 
determine whether a significantly better fit to observation 
was obtained with the latter. 

In addition to the IDL routines, cross-checking of results 
was performed using the Dipso package (part of a suite of 
tools produced for the UK Starlink system) . Excellent con- 
sistency was observed between results obtained from these 
disparate packages. 

Co-addition of spectral features has been performed 
in several cases in order to reveal details which are only 
marginally detectable at the S/N of the data. In this tech- 
nique, ~ 10 A - wide sections of spectral data are extracted, 
each centered on the wavelength of a particular line of a 
given species (e.g. the AA1548,I550 lines of CIV). The sec- 
tions are then transformed into velocity space so that each 
shows a line at the velocity of the primary interstellar cloud 
or the photosphere. Spectral sections are then co-added so 
that the strength of the (randomly distributed) noise fea- 
tures remains essentially unchanged, while absorption fea- 
tures sharing a common velocity are summed. The result is 
a spectrum with improved S/N, showing the profile of lines 
of a particular species. 

Co-addition does not improve the resolution of the data, 
and is therefore ineffective in revealing circumstellar features 
which are blended with their photospheric counterparts at 
the resolution of the data. Further, the technique is only 
effective when several of the primary lines are accompanied 
by such features. Nevertheless, co-addition has proved to be 
a useful technique in detecting weak circumstellar features 
which are clearly separated from the primary component at 
the resolution of the instrument. 



2.3 Comments on individual objects 

2. 3. 1 Stars exhibiting circumstellar features 
REJ 1738+665 

RE J 1738+ 665 is the hottest DA w hite dwarf to be detected 
by ROSAT fearstow et alJll994bT) . A photospheric velocity 
of v phot ~ 30+1 km s 1 is determined, based on absorption 
features arising from Fe, Ni and O which show no multiple 
components; interstellar lines indicate a line of sight ISM 
velocity of i>i SM ~ — 18 ± 1 kms -1 . The line of sight velocity 
of the LIC is estimated to be Uuc* ~ —3.4 kms - . 

Clear evidence is found for the presence of circum- 
stellar material in this star. Figure Q shows the C IV res- 
onance doublet in REJ 1738+665, with shifted components 
at —18.5 ± 0.5 kms -1 dominating the 30 kms -1 photo- 



REJ 1738 + 665 C IV doublet 




1547.0 1548.0 1549.0 1550.0 1551.0 

Wavelength (A) 

Figure 1. Clear multiplicity in the C IV doublet of RE 1738+665, 
with the photospheric components dominated by the blueshifted 
lines. 



spheric contribution. Shifted features with similar veloci- 
ties are observed in several other species, although in each 
case the photospheric component is dominant. The Si IV 
doublet shows non-photospheric components at —17.7 + 0.7 
kms -1 . Viewed individually, the lines of the NV dou- 
blet (AA 1238.821, 1242.804) show no evidence of compan- 
ions, but co-addition of these features suggests an extra, 
weak component at —15.2 + 2 kms -1 . The OIV doublet 
(AAI 338.612,1343.512) shows no additional features, but the 
O V line at Af 371.292 is accompanied by a weak shifted com- 
ponent at —18.7 ± 0.5 kms -1 . Curve-of-growth analysis of 
the C IV and Si IV lines indicate column densities of N (C IV) 
2 and N (Si IV) = 1.66 x I0 13 atoms 
1 (curves for the Si IV features are 



1.70 x 10 atoms cm" 



cm - with b ~ 5 kms 
shown in figure |5J. 

The velocities measured for these shifted features 
differ from Vism by less than 0.75 kms -1 . This raises the 
possibility that they may be produced by photoionisation 
of the IS M within the Stromgren s phere around the star. 
However, iTweedv fe Kwitterl il994f) also present evidence 
for the possible existence of a planetary nebula around the 
star, based on the observation of N II circumstellar features 
at optical wavelengths. The non-photospheric features of 
REJ 1738+665 may therefore be produced by ionisation 
of the ancient planetary nebula remnant surrounding this 
star. The relationship between planetary nebulas and highly 
ionised non-photospheric lines is discussed in section [3] 

Ton 021 

Other than its inclusion in general white dwarf surveys, Ton 
021 has received comparatively little attention in the liter- 
ature. We determine T eff = 69700 ± 530 K and log g = 
7.47 ± 0.05 (c.f. 6971 1 + 1030 and 7.469 ± 0.05 respectively 
from iFinlev et alJ ill 9971) ). Weighted averages of the signif- 
icant interstellar and photospheric absorption features give 
values of Wi SM = 0.86 ± 0.01 kms -1 and v phot — 37.25 ± 0.22 
kms -1 . Among the observed interstellar features are strong 
lines of CII at A1334.5323 (equivalent width 171 ± 1.6 mA) 
and AI335.7076 (equivalent width 61 ± 2.2 mA). Veloci- 
ties consistent with the weighted average Wism are observed 
for the broad interstellar lines of Nl and Sill; however, the 
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REJ 1738.665 SilV 

Log N = 13.22 



-3 - 




6 7 8 9 10 11 12 

Log(NfA) (cm"') 

Figure 2. Curves of growth for Si IV in REJ 1738+665. Separate 
curves in each plot correspond to different values of the Dopplcr 
parameter, b (indicated in kms -1 ). 



Si III A1206.5 line ap pears at 7.8 ± 0.7 km s" 1 . This is not 
a unique observation; iHolberg et ail lll999aD find significant 
differences between the velocity of this line and the mean 
value of Vis m in REJ 1032 +532 (also included in the cur- 
rent study), [rlolberg et al. I discuss this observation in some 
detail, arguing that the Si III feature is unlikely to originate 
in the LIC due to the high ionisation fraction of hydrogen 
(~ 95%) required to maintain detectable amounts of Si III, 
which has a high rate coefficient for charge exchange with 
neutral H. As noted by|HcJb^rg_e/t_alJ jl999a|) , a similar ob- 
servation is made bv IVMaT^yhid!ai^tedT7l998r) in the case 
of G191-B2B. 

Both lines of the CIV resonance doublet are accom- 
panied by shifted features. In the A1548.202 line, the non- 
photospheric component is best fit by a Gaussian with a 
velocity of w cilc = 9.5 ± 1.0 kms -1 , and an equivalent width 
of 20 ± 5 mA. This velocity agrees, within the stated er- 
ror margin, with that determined for the Si III line. The 
shifted component in the A1550.774 line is fit by a Gaussian 
at v circ = 5.5 ± 2.2 kms -1 , with an equivalent width of 9.4 ± 
2.5 mA. The estimated CIV column density contributing to 
these shifted features is N (CIV) = 3.98 x 10 12 - 1.26 x 10 13 
atoms cm -2 based on a curve-of-growth analysis. In both 
cases, the photospheric components are found at velocities 
consistent with the average value for v vllot . The presence of a 
non-photospheric CIV feature is most clearly demonstrated 
in the co-added lines of the resonance doublet, as shown in 
figure |3 

Viewed individually, the lines of the Si IV resonance 



Ton 21 co — added C IV features Ton 21 co— added Si IV features 




-40 20 60 100 -40 20 60 100 

Velocity (km/s) Velocity (km/s) 



Figure 3. The co-added C IV and Si IV resonance lines of Ton 021. 
Clear evidence for multiplicity is observed in the C IV features, 
while the presence of a non-photospheric component to the Si IV 
features is less compelling. Dotted line: best fit dual-Gaussian to 
the data. 

doublet appear slightly asymmetrical but show no clear 
multiplicity at the resolution of this data. When co-added, 
this asymmetry is more noticeable (as shown again in 
fig CUi an d we find the best fit dual-gaussian to be one 
with a primary component at ~ 37 ± 1.2 kms -1 , in 
agreement with v phot , and a non-photospheric component 
at ~ 6.5 ± 12.3 kms -1 . Although this velocity is close 
to the value of u circ found in the CIV doublet, there is 
considerable uncertainty in the measurement, and data of 
higher S/N will be required to confirm the detection of 
non-photospheric components to the Si IV doublet of Ton 
021. 

REJ 0558-373 

The CIV AA1548. 202, 1550. 774 lines in this star are moder- 
ately asymmetrical. An F-test suggests that a dual Gaussian 
fit is preferred over a single line, above the 98% confidence 
interval, for each feature. No corresponding asymmetries are 
observed elsewhere in REJ 0558-373. 

For the A1548.202 line, the individual components are 
found at 7.0 ± 1.0 and 26.8 ± 1.2 kms -1 (with equivalent 
widths of 79 and 68 mA respectively). Corresponding values 
for the A1550.774 line are 7.9 ± 1.2 and 21.8 ± 1.0 kms -1 
(13 mA and 150 mA). Fits to each line are illustrated in 
figure 2] For comparison, u ph ot = 22.7 ± 2.8 kms -1 , and 
Vism = 11 -6 ± 1.4 kms -1 , and hence the longer wavelength 
component of each C IV line is in reasonable agreement with 
the photospheric value. 

fgrav is estimated as ~ 20 kms -1 for this star, which 
is greater than the ~ 10 kms - difference between the 
photospheric and shifted C IV components. It is therefore 
possible that the blueshifted, non-photospheric CIV fea- 
tures could be formed by material within the potential well, 
rather than the weakly shifted outer regions. However, the 
non-photospheric components lie close to the velocity of the 
ISM, raising instead the possibility that the star is ionising 
material in its local interstellar environment. In either case, 
the absence of corresponding features in other strong lines 
is puzzling. 
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Figure 4. Lines of the C IV doublet of REJ 0558-373, in velocity 
space, with Gaussian fits overlaid (histogram = STIS data; dashed 
curves: component Gaussians; solid curve: compound fit). 



WD 2218 + 706 co-added C IV doublet 
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Figure 5. Histogram: co-added lines of the CIV doublet in WD 
2218+706 showing the photosphcric component at 38.8 kms , 
accompanied by a redshifted feature at -16.3 kms -1 . The data 
are accompanied by the individual and summed components of 
the best-fit dual Gaussian. 



WD 2218+706 

This object is unusual and important in several respects. 
Lines of the CIV and Si IV doublets are clearly multiple, 
dominated by a photospheric contribution, but with accom- 
panying components of comparable equivalent width at a 
velocity of -16.3 ± 0.7 kms -1 . These features are therefore 
redshifted with respect to the photospheric velocity (v phot = - 
38.7 ± 0.2 kms -1 ), possibly representing the infall of ma- 
terial onto the white dwarf. Gravitational redshift therefore 
provides no viable explanation for these lines. 

WD 2218+706 is surround ed by an old planetary neb- 
ula, D eHt5, and is discussed bv iNapiwotzki fe Schonbernerl 
lll99 5 j). In a st u dy o f planetary nebula dynam- 
ics. Ipgani fc Sokerl (Il998l) show that in regions where 
the ISM is reasonably dense (such as the galactic plane), 
Rayleigh- Taylor instabilities can develop in the outer 
regions of planetary nebulae, leading to fragmentation of 
the halo, and allowing the surrounding ISM to pass into 
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Figure 6. Curves of growth for the redshifted C IV features in 
WD 2218+706. Separate curves in each plot correspond to differ- 
ent values of the Doppler parameter, b (indicated in kms -1 ) 



the inner regions of the nebula where photoionisation can 
occur. Although WD 2218+706 is out of the galactic plane 
(fen = 11.6°), and therefore lies in a region where the mean 
ISM d ensity may be expected to be relatively low, iKunl 
Jl998f) describes the morphology of a nearby giant molecu- 
lar cloud complex consisting of a large number of distinct 
regions previously identified in independent surveys; several 
are found close to WD 2218+706, and two are of partic- 
ular interest (Lynds 1217 and Lynds 1219). The central 
portions of these clouds have galactic coordinates within 
less than 0.5° of this star, and their distance limits (from 
380 to 450 pc ) encompass the distance t o WD 2218+706 
(440 pc, from INapiwotzki fc Schonbernerj) . This raises the 
possibility that the star may lie in an area where the ISM 
is particularly dense, allowing instability and inflow to take 
place (a curve of growth analysis for the non-photospheric 
features in WD 2218+706 suggests column densities of N 
(CIV) = 4.17 x 10 13 atoms cm -2 and N (Si IV) = 4.07 x 10 13 
atoms cm -2 , each with a Doppler parameter of 6 kms -1 ). 
However, alternative explanations, such as the presence of a 
hidden companion, are also deserving of investigation, and 
this work is currently in progress. 

During the course of the WD 2218+706 study, evi- 
dence was found for t he existence of trace amounts of He 
in the STIS spectrum fearstow et al]|200ll) . with the Hell 
A1640.5050 line observed in the STIS spectr um close to the 
estimated photospheric velocity. As noted bv lBarstow et all 
the A1640.5050 line has an n=2 lower level, which should not 
be populated in collisionless interstellar material, while any 
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helium in the surrounding planetary nebula would be ex- 
pected to be found in emission. Hence a photospheric origin 
appears to be the most satisfactory of these three possible 
sources. The surface gravity of WD 2218+706 (log g w 7.00) 
is low for an isolated white dwarf, and may be explained 
in terms of close-binary evolution, in which the progenitor 
star fills its Roche lobe and loses mass to a companion. This 
loss of material prevents helium ignition from taking place. 
Instead, the star, consisting of a He core surrounded by a H- 
rich envelope (which still supports nuclear reactions at the 
base), contracts slow ly towards the low-mass, He-core whit e 
dwarf configuration jDriebe et alJll998t lrfapiwotzkilll999l) . 
As the case of Feige 24 illustrates, the presence of a binary 
companion can contribute to the appearance of circumstellar 
lines. However, there is no direct evidence for the existence 
of a binary companion, and the range of possible masses for 
this star does not preclude AGB evolution, leaving open the 
possibility that WD 2218+706 is the product of single-star 
evolution. 

Evidence exists for at least one other H-rich white 
dwarf star exhibiting the Hell A1640 line, in the DAB 
HS 0209+0832 (T c ff ~ 35,000 K; l og g ~ 7.8). This 
star is discussed bv IWolff et alJ l)2000F) . who suggest that 
significant quantities of He are present in the atmosphere, 
despite the short diffusion timescales for He, as a result 
of ongoing accretion of matter from an interstellar cloud. 
Supporting evidence for this e xplanation can be found in 
the work bv lHeber et aJ l)l997f) . who observe variability in 
the strength of the He A1640 line, possibly as a result of 
the passage of HS 0209 +0832 through an inho mogeneous 
medium. Alternatively, lUnelaub fc Buesl l|2000l) find that 
DAO stars, in which mass loss in the form of a stellar 
wind prevents He from sinking out of the atmosphere, can 
transform into DA stars when the phase of wind-driven 
mass loss ends. This transition is found to occur near 
log q ~ 7.0 for a star with T cfr ~ 60,000 K (see fig.6, 
llJnglaub fc Buesl ^2000T l ^1 . raising the possibility that WD 
2218+706 may be such a transitional object. 

Feige 24 

Feige 24 is a white dwarf+red dwarf binary sys- 
tem, a nd is the subject of s everal detailed stud- 
ies (e .g. iDupree fc Raymond! Il982l IVennes fc Thorstensenl 
119941) . Two STIS data-sets were available for this star, ac- 
quired on November 29 th 1997 (binary phase 0.73-0.75) and 
January 4*^ 1998 (binary phase 0.23-0.25), repres enting the 
orbital quadrature points. fVennes fc Thorstensenl estimate a 
systemic velocity of 62.0 ± 1.4 km s - . For the current study, 
the systemic velocity has been estimated by taking the mean 
of the photospheric values obtained from each data set (31.6 
and 129.1 kms -1 ), resulting in an estimated systemic ve- 
locity of 80.3 ± 0.5 kms -1 . Uism is estimated at 8.2 ±0.1 
kms -1 . 

Feige 24 is known to exhibit multiple components, in the 
lines of the C IV doublet only. In this work, the dominant 
components match the photospheric velocity of each data- 
set, and secondary features are observed to remain at 7.8 ± 
0.2 kms -1 , irrespective of the orb ital phase. This statio nary 
component has been discussed bv lDupree fc Raymond! who 
suggest that the most likely source is a Stromgren sphere ex- 
cited by the white dwarf, and measure column densities of 
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Figure 7. The C IV 1550 A feature of Feige 24 as observed at the 
two quadrature points. The stationary circumstellar component 
is clearly visible, as is the orbital velocity of the photospheric 
component. 

N(CIV) = 3.86(±1.51) x 10 13 a toms cm -2 using curves-of- 
growth. IVennes fc Thorstensenl investigated the possibility 
that the material responsible resides in the photosphere, in 
a circumstellar shell, or in a wind from the red dwarf com- 
panion. The current study shows no shifted components in 
any of the other resonance lines (e.g. Si IV or N V). A third, 
very weak feature on the red side of each photospheric line 
in the CIV doublet, most obviously at 1550 A, is identified 
as the FeV A 1550.907 line. 

The gravitational redshift estimate made during this 
study, D„, v sa 17 kms -1 , is som ewhat higher than that 
derived bv IVennes fc Thorstensenl (u grav ~ 9 ± 2 kms -1 ), 
but is still too low to explain the secondary CIV compo- 
nents. However, their velocities agree, within error, with 
that of the ISM, and hence a link between the star and its 
immediate surroundings (beyond a ny circumstellar shell) 
cannot be discounted. Alternatively. [Vennes fc Thorstensenl 
estimate that a C IV column density from N (C IV) = 
7.94 x 10 11 - 7.94 x 10 13 atoms cm (corresponding 
to equivalent widths of between 4-400 mA for the non- 
photospheric component of the A 1550.774 line, assuming a 
linear curve-of-growth) , would be consistent with mass loss 
from the red dwarf companion. Although insufficient data 
are available to derive an unambiguous C IV column density 
in this study, the estimated value of N (CIV) = 1.48 x 10 13 
at oms cm -2 , is within the range of possible values obtained 
bv lVennes fc Thorstensenl and the equivalent width of the 
A 1550.774 line (24 mA) is also consistent with the large 
range allowed by earlier estimates. 

G191-B2B 

The STIS E140M data for this star show the CIV 
resonance doublets are accompanied by strong non- 
photos pheric blue shifted com ponents, confirming the re- 
sults of lBruhweiler et alJ (I1999T) . The wavelength and equiv- 
alent width of these features have been estimated by fit- 
ting multiple Gaussians to the data, as shown in figure |H| 
The component at photospheric velocity is of greater equiv- 
alent width, and the shifted element s are observed atve- 
locities of 7.6 ±0.2 kms -1 . Recently. iHolberg et all (I2002T) 
also detected weak blue shifted components to the Si IV 
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Figure 8. The C IV doublet of G191-B2B, in velocity space, with compound Gaussian fits (solid lines), and the individual Gaussian 
components (dashed lines). 



AA1393. 755, 1402. 777 doublet at the same velocity as the 
non-photospheric CIV features, using 22 co-added STIS 
E140H spectra. 

The calculated figure of w grav » 15.4 kms -1 is com- 
parable with the velocity difference between photospheric 
and shifted high ionisation features, suggesting that the non- 
photospheric material probably resides outside the limit of 
the potential well. The velocity of the highly ionised non- 
photospheric features is substantially different to the value 
of the interstellar features («ism = 16 ± 1 kms -1 as de- 
termined in this work), and at first sight, photoionisation 
of the cloud responsible for the primary ISM features does 
not appear to provide a viable explanation. However, the 
value quoted for Vi SM is based on analysis of an E140M 
(medium re solution) STIS dat a-set, with a resolving power 
of ~ 35000. ISahu et alj (Il999l) describe observations made 
with the E140H grating (resolving power ~ 110000), and 
clearly show two distinct interstellar components, with ve- 
locities of ~ 8.6 kms and ~ 19.3 kms -1 , the latter com- 
ponent having a velocity close to the predicted value of « L ic 
(estimated at 20.58 kms -1 in this study). Clearly, the highly 
ionised non-photospheric components have a velocity which 
is very close to the 8.6 kms - interstellar cloud. 

A curve-of-growth analysis was performed for the CIV 
features in G191-B2B. As in the case of Feige 24, the 
availability of only two datum points prevents any rigorous 
constraints from being placed on the implied CIV column 
density, though the value of N (CIV) = 2. 40 x 10 13 atoms 
cm" 2 is not dissimilar from the results of IVennes fc Land 
teOOlf) . who use synthetic modeling techniques to estimate 
a value of N (CIV) = 6.31 x 10 13 atoms cm" 2 . The Doppler 
parameter suggested by the current analysis, b = 10 kms -1 , 



is significant ly higher than the value of b = 5.2 kms - 
presented bv lVennes fc Lan3 this discrepancy may also be 
explained by the lack of available data in the current work. 

REJ 0457-281 

The exceptionally low HI column density to this star 
(1.3 x 10 17 atom cm 2 ), was revealed in the discovery pa- 
per bv iBarstow et al.l Jl994al) . Along with G191-B2B, this 
white dwarf was the first t o have phosphorou s and sulphur 
identified in its spectrum llVennes et al.llT996l) . Later, HBS 
showed that the photospheric Si IV and C IV resonance lines 
of REJ 0457-281 are accompanied by blue-shifted features 
(figure|5J. 

Few interstellar and photospheric lines are identifiable, 
making precise velocity measurements difficult. The ISM ve- 
locity estimate of HBS is confirmed, but somewhat higher 
photospheric velocities are derived. Co-addition of the CIV 
doublet lines clearly reveals two velocity components: one 
at 22.5 ±1.39 kms -1 , and the photospheric component at 
81.26 ± 2.65 kms -1 . A curve of growth analysis for the non- 
photospheric components suggests N (CIV) ~ 1.82 x 10 14 
atoms cm -2 with a Doppler parameter b « 4 kms -1 . For the 
co-added SilV doublet, corresponding velocities are 19.08 
± 4.31 kms -1 and 80.65 ± 1.38 kms -1 . Although multi- 
ple velocity components are not obvious in the NV dou- 
blet, there is, nevertheless, evidence to suggest that they 
are present (figure|UJ. The A 1238.8210 N V line is a narrow, 
well defined feature at 76.33 ± 4.78 kms -1 , accompanied 
by a weaker blueshifted feature at 16.55 ± 4.61 kms - . The 
main A 1242.804 line has a similar velocity (76.91 ± 4.38 
kms -1 ), but shows only tentative evidence for a blueshifted 
component. 
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Figure 9. From top to bottom, the co-added lines of the N V, 
Si IV and C IV doublets in REJ 0457-281. Lower arrow indicates 
the approximate photospheric velocity, upper arrow shows the 
position of the shifted components, present in each species. 



From these data, a weighted average is computed for 
the photospheric and blueshifted velocity components, 
suggesting v phat = 76.91 ± 0.83, (c.f. 69.60 ± 1.97 from 
HBS), and « circ = 21.76 ± 1.27. Thus, the estimated velocity 
shift of the blueshifted features relative to the photospheric 
components agrees, within the stated error, with the 53 
kms -1 value of HBS. 

REJ 2156 - 546 

iBarstow et alJ il997T) determined limits to the heavy ele- 
ment abundance in REJ 2156-546, and described this ob- 
ject as being similar to HZ 43 in having a reasonably pure 
H atmosphere. The STIS spectrum of REJ 2156-546 shows 
clear interstellar lines, indicating Vjsm = -8.39 ±0.17 kms - . 
These new, high resolution data also appear to show features 
due to photospheric material. The lines are weak, and unam- 
biguous identifications are limited to the strong resonance 
doublets of Si IV and CIV, although there may be features 
from NV and NiV at the detection limit. To obtain a re- 
liable value for v phot , the Si IV lines (AA1393. 755, 1402. 777) 
were co-added in velocity space, producing a clear feature at 
-17.79 ±1.33 kms -1 , with an equivalent width of 11 mA. Co- 
addition of the C IV doublet appears to reveal two features 
(figure Hill. The first, at -20.71 ±0.80 kms -1 , is weak (4.7 
mA), and very close to our value for v P h Q t- The dominant 
second component lies at -1.65 ±0.76 kms -1 , and has an 
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Figure 10. Co-added Si IV and C IV features in REJ 2156-546, 
with Gaussian fits overlaid. 



equivalent width of 20 mA. The weighted average of puta- 
tive photospheric features produces a value of v phot = -19.94 
±0.68 kms -1 . 

The proposed CIV feature lying at a similar velocity 
to the single Si IV line is admittedly weak, and spectra 
of improved S/N will be required before these results 
can be regarded as incontrovertible. Nevertheless, if it is 
assumed that the object is devoid of any non-photospheric 
features, whether blue- or redshifted, then the relatively 
large difference in velocity between the Si IV line and 
the dominant CIV feature (approximately 16 kms -1 ) is 
somewhat difficult to explain. Clearly, this is an object 
deserving of further attention. 

REJ 1611 - 085 

iHoIberg et alJ il997t) found the amount of Si in the spec- 
trum of REJ 1614-085 to be an order of magnitude under- 
abundant compared to the predictions of radiative levitation 
calculations, while N appears to be three orders of mag- 
nitude over-abundant. Two velocity components were ob- 
served in the line of sight ISM, but most significant for the 
current work is the result that the lines of the C IV and Si IV 
doublets exhibit weak blueshifted features, as illustrated in 
figure HT1 

Results from this study suggest that the primary ISM 
component lies at a velocity of -29.56 ± 0.33 kms -1 , and 
the secondary component at ±48.64 ± 1.15 kms -1 . These 
values compare reasonably well with those of lHolberg et all 
who find velocities of -27.05 ± 1.5, and ±47.40 ±1.50 
kms -1 respectively. Similar agreement with the earlier 
study is also found when considering the photospheric and 
blueshifted features. The photospheric velocity is found to 
be fphot = -37 .31 ±0.40 kms -1 , in excellent agreement with 
IHoIberg et all As recorded previously, no clear evidence 
exists for a secondary component to the photospheric 
NV lines, although it should be noted that the centroid 
positions of Gaussian fits to these lines are found to differ 
by approximately 5 kms . In the case of the SilV doublet, 
the secondary components are shifted by -40 kms -1 relative 
to the primary features, as determined by IHoIberg et alJ 



1997); for the CIV doublet, this figure is -29 kms , com 



quoted by IHoIberg' et alJ 



pared to the value of -25 kms 
a 1997T) . This apparent discrepancy is most likely to be a 
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Figure 11. Co-added lines of the C IV (left) and Si IV (right) 
doublets in REJ 1614-085, in velocity space. The bl ueshifted com- 
ponents previously noted by iHolberg et all Il997l) are indicated 
by arrows. 

result of the different positions chosen for line demarcation 
in the two studies. A curve of growth analysis performed 
on the shifted CIV features suggests N (CIV) « 3.16 x 10 13 
atoms cm" 2 , and the Doppler parameter b &s 2 kms , 
though with only two data points, these values are partic- 
ularly poorly constrained, and must not be over interpreted. 

GD 659 

Both IUE and STIS data are available for this star, although 
the STIS spectrum is of limited coverage (1160 - 1357 A). 
The ISM velocity determined from the IUE data (wism = 
12.33 ± 1.52 kms" 1 ) agrees with that of HBS, while the 
IMS-based photospheric velocity appears to be somewhat 
lower than previously quoted (v ptlo t = 33.51 ± 1.03 kms - , 
c.f. 40.31 ± 1.83 kms" 1 from HBS). 1 

Since these measurements were made with identical 
data sets, this discrepancy must be ascribed to differences 
in choices of continuum levels and line boundaries used dur- 
ing the measurement process. However, the available STIS 
data also point to a lower photospheric velocity, with v ptlat 
= 34.28 ±0.17 kms -1 , in agreement with the IUE estimate 
made in this st udy, and close to the value of 33.58 kms" 1 
determined by IHolberg et ail feOOd) using STIS data. The 
STIS ISM velocity is also lower than the IUE value, at 
Wism = 9-77 ± 0.22 kms" 1 . These STIS velocities, based on 
higher resolution data with better S/N, are adopted for GD 
659 in tabled Although the resonance doublets of N V, Si IV 
and C IV are clearly visible in the IUE data, the resolution 
is insufficient to observe well defined Gaussian profiles, and 
thus the sensitivity to any non-photospheric components is 
low. However, the profiles appear to be narrow, ruling out 
any obvious multiplicity in these lines. STIS data show the 
lines of the NV doublet as narrow and symmetrical, effec- 
tively ruling out the existence of non-photospheric N V com- 
ponents. 

1 These features were originally attributed to circumstellar ma- 
terial, and are listed accordingly by HBS, based on the 76 kms 
velocity difference between the UV feat u res an d the -37 kms 
Balmer line radial velocity of IWeenerl <1974D . The blueshifted 
IWegnerl value is most likely the result of a typographical error. 
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Figure 12. Left: the A 1548.202 (top) and A 1550.774 (bottom) 
C IV features in GD 659. Dashed vertical lines indicate the posi- 
tion of a possible secondary component near kms -1 , and the 
primary photospheric line. Flux levels have been adjusted by a 
constant to allow plotting of features in the same panel. Right: 
co-addition of the C IV features, with a best-fit double Gaussian 
profile overlaid. 

One interesting feature, clearly visible in the CIV 
A 1548.202 line, though present also in the A 1550.774 
line, is a weak feature near kms -1 (figure [T^ll . Fitting a 
double Gaussian to the co-added CIV lines, velocities of 
36.74 ± 2.56 kms" 1 (primary), and -2.97 ± 3.00 kms" 1 
(secondary) are obtained. The secondary component is 
weak (with an equivalent width of 6 mA compared to 36 
mA for the photospheric component), and is of comparable 
strength to scatter in the adjacent continuum regions. 
Although an F-test indicates that a dual Gaussian fit is 
preferred over a single component at the 94% confidence 
interval, the similarity between this feature and the natural 
scatter in the data suggests that this is simply noise. 
However, until high resolution STIS data for this region 
can rule out the existence of shifted components in the C IV 
lines, we tentatively include GD 659 among the stars with 
possible circumstellar features. 



2.3.2 Stars with no clear circumstellar features at the 
resolution of current data 

PG 0948+534 

PG 0948+534 is the hottest star in this sample. Strong ISM 
lines are observed, many of which are saturated, and at least 
two velocity components are present. However, Nl and Sill 
features are most accurately described by three components, 
with self consistent velocities: the primary component lies at 
-0.26 ± 1.26 km s" 1 , with the secondary and tertiary compo- 
nents at -22.8 ± 1.2 and 22.1 ± 2.4 kms" 1 . These are most 
clearly revealed by co-addition of the Nl and Sill features, 
as shown in figure ED 

The photospheric lines, including the resonance dou- 
blets of CIV and Si IV, exhibit narrow, symmetrical profiles 
and are apparently devoid of any shifted components, 
defining v phot — -14.25 ± 0.22 kms" 1 . A remarkably 
strong, multi-component CII 1335.7076 A feature is ob- 
served, also shown in figure 1131 The velocity components 
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Figure 13. Upper left: co-added interstellar Si II lines in 
REJ 0948+534, with three-component Gaussian fit. Lower 
left: analogous plot for N I. Right: interstellar C II lines at 
AA1334.5323,1335.7076. 



Figure 15. Left: lines of the C IV doublet in REJ 0623-371, 
in velocity space, demonstrating the difference in velocity of the 
two lines (A1548 = solid histogram, A1550 = dashed histogram). 
Right: dual Gaussian fit to the 1550 A C IV line (dashed curves 
= component Gaussians, solid curve = summed Gaussian. 
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Figure 14. Observed C IV doublet in REJ 2214-492 (black) 
compared to a synthetic spectrum with N(C)/ = 10" 6 N(H), 
smoothed to the resolution of IUE. 

match those of other ISM lines, although the -23 kms" 1 
feature is very weak, manifesting itself as a broadening 
on the blue side of the line. Excited Sill transitions 
such as AA1265.002 1309.276,1533.431 are not observed 
llHolberg et alJ <ll995l) used the presence of these lines in 
the white dwarf CD -38° 10980 to infer the existence of a 
circumstellar cloud around the star). No evidence of highly 
ionised non-photospheric material is found in PG 0948+534. 

REJ 2214-492 

Weighted average line velocities indicate vi SM = -1.72 ± 0.51 
kms" 1 , and v phot = 33.49 ± 0.45 km s - . These values com- 
pare well with those of HBS, who find Vism = -0.71 ± 0.88 
and Vphot = 33.91 ± 0.47 kms" 1 respectively. However, a 
significant difference exists between the velocity of each line 
in the CIV doublet, with A 1548.202 at 30.5 ± 2.1 kms" 1 , 
and A 1550.774 at 40.4 ± 2.7 kms" 1 , if each of the two lines 
is assumed to be made up of only one absorption feature. 

Visual inspection of the CIV doublet reveals slight 
asymmetry, particularly in the A 1548.202 line, where a 
dual fit was found to be superior to the single Gaussian 



at the 99.9% confidence level, with velocity components at 
5.37 and 38.36 kms -1 , and equivalent widths of 31.1 and 
99.3 mA respectively. A dual Gaussian fit to the A 1550.774 
line produced a less obvious improvement (at the 90% 
confidence level) with components at 36.2 and 80.6 kms - , 
and equivalent widths of 137.6 and 14.8 mA respectively. 
The primary Gaussian components of the doublet thus lie 
at velocities more consistent with each other and the overall 
photospheric value. The status of putative non-photospheric 
contributions in the CIV doublet is less certain; although 
the 5 kms -1 feature at 1548 A appears to provide a 
good match to observation, the lack of a corresponding 
feature at 1550 A prevents confirmation of its reality. No 
evidence was found for multiplicity in other photospheric 
lines. Line profiles were compared with those from a model 
spectrum, produced using the TLUSTY and SYNSPEC 
codes, and a dopting the heavy el ement abundances de- 
termined by iBarstow et"afl fcOOlfl (with N(C)/N(H) = 
1.0 xlO -6 ). After smoothing the model to the resolution 
of IUE, no significant differences were apparent in the 
shapes of model and observed CIV lines (figure [TH . casting 
further doubt on the presence of shifted features in this star. 

REJ 0623-371 

The photospheric velocity determined here (v phot — 41.15 
±0.56 kms" 1 ) agrees with that of HBS, though the cur- 
rent value for Vism is somewhat lower than given by HBS 
(16.40 ±0.70 kms" 1 c.f. 19.48 ±0.85 kms" 1 respectively); 
the principal source of this difference lies in a more precise 
determination of the C II line velocity. HBS estimate the ve- 
locity of this line to be 14.07 ±3.24 kms" 1 , compared to the 
new value of 13.33 ±1.14 kms -1 . 

No compelling evidence exists for the presence of 
shifted components in the spectrum of REJ 0623-371, but 
as in the case of REJ 2214-492, a significant difference is ob- 
served in the velocity of the lines in the CIV doublet (1548 
A = 38.6 ±2.3 kms" 1 , 1550 A = 47.6 ±2.2 kms -1 ). In 
contrast, the lines of the NV and Si IV resonance doublets, 
which are of comparable equivalent width, agree within the 
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estimated error. Determining the reality of any features in 
the doublet is complicated by considerable absorption in 
the continuum of this extremely metal rich DA, though by 
restricting Gaussian fits to the region below this structure, 
useful comparisons between the level of agreement found 
with single and double Gaussian profiles, may be obtained. 
Using this method, a dual Gaussian fit is preferred to a 
single feature only at the 88.8% level for the 1548 A line. 
However, results for the 1550 A line are less ambiguous, 
suggesting a dual fit at 98.9%. The resulting Gaussians have 
velocities of 45.1 and 77.3 kms" 1 , and equivalent widths 
of 128.5 and 10.8 mA respectively (figure ^J. Thus, while 
the quality of the available data is insufficient to prove the 
existence of circumstellar features in the star, the results of 
this analysis provide some justification for proposing repeat 
observations at a higher signal-to-noise ratio (S/N) and 
spectral resolution. 

REJ 2334-471 

Values measured for Vism and v pilat are in good agreement 
with those obtained by HBS. The relatively poor quality of 
the IUE data precludes unambiguous identification of any 
non-photospheric features, particularly in the case of the 
CIV lines. There is no evidence of multiplicity in the NV 
doublet, although a dual Gaussian fit to the 1242 A line, 
with components at 19.7 and 43.51 kms" 1 , produces a fit 
which is preferred over a single feature at the 93% confidence 
level. 

It is therefore intriguing that each line in the Si IV 
doublet (AA1393.755, 1402.777) is fitted reasonably well 
(above the 95% confidence interval when compared to 
a single feature) by double Gaussian profiles. Each line 
can be described by a double Gaussian with Vi = 34.00 
±0.82 kms" 1 , EWi = 48.85 ±0.95 mA, and V 2 = 54.64 
±1.58 kms" 1 , EW 2 = 22.74 ±3.25 mA. Neither of the 
Gaussian velocities are in agreement with the average 
photospheric value, although a considerable spread in 
individual photospheric velocity measurements is observed, 
so that the discrepancy cannot be used to infer the absence 
of such features. Spectra with improved S/N are required to 
confirm or disprove the existence of circumstellar features 
in this star. 

GD 246 

IUE and STIS data were available for this star. Photo- 
spheric and ISM line velocities measured from the IUE data 
are in good agreement with those of HBS. The photospheric 
velocity estimated from STIS data differs from the IUE 
value by 1 kms - , although this is within the IUE error 
bounds. The STIS ISM velocity (-5.78 ±0.12 kms" 1 ) is 
marginally outside the IUE estimate (-7.87 ±1.00 kms" 1 ). 
IUE and STIS data clearly show the C IV doublet lines as 
singular, and at the photospheric velocity. STIS data shows 
the Si IV A 1393.755 line as being devoid of any secondary 
components, while a sharp feature (only one bin, or 0.02 A 
in width) is observed on the redward edge of the A 1402.777 
line. Although IUE data also hint at a broadening on this 
side of the line, unconvincing dual Gaussian fits, and the 
extreme narrowness of the extra feature in STIS data, 
suggest that this is due to noise, and thus GD 246 shows 
no clear evidence of circumstellar material. 



PG 1123+189 

In the photometric study of hot white dwarfs bv lGreen et alJ 
i2000T) . this object is one of those listed as having a sig- 
nificant IR excess, suggesting the possibility of a low mass 
companion to the star. 

The STIS spectrum for this object is limited in coverage 
(1163 - 1361 A), and has relatively low S/N. However, 
many interstellar lines are visible in the spectrum, and a 
value of Vism = -0.67 ± 0.06 kms" 1 is obtained. Although 
photospheric features are difficult to distinguish in the 
data, by co-adding the NV lines at 1238 and 1242 A with 
those of NiV between 1250 and 1336 A in velocity space, 
a single absorption feature is clearly visible, suggesting a 
value of Uphot = 12.55 ± 0.53 kms" 1 . The quality of these 
data is insufficient to confirm or rule out the presence of 
non-photospheric features with confidence. 

HZ 43 

A well studied object, HZ 43 is a member of the group of 
white dwarfs which can be adequately modeled with an 
atmosphere devoid of any heavy elements. Several ISM 
lines are observed, leading to an estimate of Ui S m which 
agrees with that of HBS. Co-addition of the spectrum at 
the wavelengths of the major N, C, Ni and Si lines fails to 
reveal any photospheric features. Similarly, co-addition, in 
velocity space, at the wavelengths of the excited Si tran- 
sitions (AA1264.738, 1265.0020, 1309.2758 and 1533.4312) 
also shows no new features. 

RE J 1032+532 
Thi s object is t h e sub j ect of a com p rehensi ve study 
by IHolberg et all (I1999T) : iHolberg et all (Il999al) . In the 
current work, the measured value of the primary ISM 
features, Vism = 0-84 ±0.21 kms" 1 , agrees, within error, 
with that of IHolberg et alJ dl999aT) . A previously noted 
secondary component to the Sill AA1193.2897, 1260.4221 
and 1526.7065 lines is found to have a velocity of -30.43 
±1 .39 kms" 1 , al so in agreement with the value quoted 
by IHolberg et all A value of u phot = 38.16 ±0.40 kms" 1 
is determined for the photospheric velocity. The excited 
Sill lin es found around som e stars possessing circumstellar 
clouds (IHolberg et alJll995l) are absent, and in none of the 
photospheric lines is any compelling evidence found for the 
existence of secondary components. 

PG 1057+719 

This object (alternative ID REJ 1100±713) is also included 
in the photometric study of white dwarfs by iGreen et alJ 
l)2000Fl . with no significant IR excess being detected. It be- 
longs to the low opacity metal poor class which inc ludes 
the majority of DA white dwarfs. THolberg et alJ lll997l) pre- 
sented a study of this star and REJ 1614-085 (see below). 
Their results revealed no signs of circumstell ar features. 

The current work confirms the results of lHolberg et all 
revealing no shifted features in the GHRS data. Co-addition 
of the ISM lines results in «i SM = -2.89 ±0.69 kms" 1 . As 
expected for a low EUV opacity object, no significant 
photospheric lines are observed. To detect any signs of 
photospheric features, a series of 10 A - wide sections 
were extracted from the data, each centred on the rest 
wavelength of one of the lines of the NV, CIV and Si IV 
resonance doublets. The sections were then transformed 
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into velocity space, and co-added. The presence of barely- 
detectable quantities of N, C and Si might then be expected 
to produce a noticeable reduction in continuum level 
around the photospheric velocity. The co-added data does 
indeed reveal a feature, with a velocity of 75.35 ± 2.59 
kms -1 , consistent with the value of 

^phot = 76.1 ±3 

kms -1 deter mined from B almer line fitting. However, as 
indicated by iHolberg et al.l weak individual features found 
near the expected positions of these lines show a consid- 
erable spread in velocity, casting doubt on their authenticity. 

GD 394 

GD 394 is photometrically variable in the EUV, though no 
signs of spectroscopic variation have been detected. In con- 
trast to REJ 1614-085, GD 394 has an extr eme overabun- 
danc e of Si compared with m odel predictions iHolberg et alJ 
Il997l) . iDupuis et alJ i2000F) note that this extreme Si abun- 
dance, and the observed EUV variability, give a unique sta- 
tus to GD 394. IDupuis et alJ present spectroscopic and tim- 
ing analyses of GD 394, which suggest the presence of a 
large EUV-dark spot on the surface of the star, sharing the 
stellar rotation period of 1.150 days. Episodic accretion is 
proposed as the source of this spot, with a magnetic field 
directing material onto the magnetic poles. No evidence ex- 
ists for the presence of a magnetic field in GD 394, though 
only upper limits can currently be placed on the strength of 
any such field. GD 394 appears to be an isolated star, and 
hence no obvious candidate exists for the source of accreted 
material, other than the immediate stellar neighbourhood. 

Early results suggested that the velocity of Si III and 
Si IV lines differed considerably from the established radial 
velocity, a nd that these lines were the refore of a circumstel- 
lar origin feruhweiler fc Kondoll983ft . It was also suggested 
that the absence of any observable Si II features in IUE data, 
which models predicted would be present, represented fur- 
ther evidence for t he non-photospheric nature of the heavy 
elements. However. lBarstow et alJ dl99rJ) demonstrated that 
the previous radial velocity measurement, obtained from the 
Balmer lines, was in error, and a revised value more consis- 
tent with the velocities found for the Si features was ob- 
taine d. By using the lates t NLTE models available at the 
time. lBarstow et alJ lll996l) showed that the predicted abun- 
dance of photospheric Si would yield line strengths within 
the noise of the IUE spectrum, and hence the absence of 
Sill features did n ot require a non-photospheric solution. 
IChav er et alJ i2000fl report the first firm detection of heavy 
elements other than Si in the spectrum of GD 394, with spec- 
tra from the Far Ultraviolet Spectroscopic Explorer (FUSE) 
showing lines of Felll and PV in the photosphere, as well 
as a large number of Si III and Si IV lines. 

Values of v 1SM = -7.28 ± 1.42 and v phot = 28.75 ± 0.91 
are obtained in this work for the velocity of I SM and photo- 
spheri c features respectively. As reported bv lHolberg et all 
il997Tl . the inventory of photospheric lines is dominated 
by S i HI, with addi tional contributions from the Si IV dou- 
blet. ^ojbe^^^^L| also observe photospheric AlIII lines at 
AA1854.7159 and 1862.7900. However, CIV and N V are con- 
spicuous only by their absence; co-addition of the data for 
these wavelengths reveals no trace of features. 

No record exists, in recent literature, of circumstellar 
features in the spectrum. The IUE data reveal no sign 
of CIV or NV, and GHRS data are available only for the 



ranges 1290 - 1325 A and 1383 - 1419 A, which do not 
include these ions. However, many Si lines are visible in 
the higher resolution GHRS data, and neither individual 
features or the co-added profiles of all visible Si lines 
show any form of asymmetry or other qualities which may 
indicate the presence of circumstellar features. GD 394 is 
therefore one of the few stars in this survey for which it can 
be conclusively stated that no shifted features exist, at the 
resolution of currently available data. 

GD 153 

Another example of a star which may be modeled with a 
pure-H atmosphere, GD 153 is a frequently observed stan- 
dard star. No obvious photospheric features are observed 
in the IUE spectrum of this star, although several ISM 
lines are recorded, indicating a value of vism = -8.42 ± 2.68 
kms -1 . This velocity agrees with that obtained by HBS. 

EG 102 

IHolberg et"afl (Il997f) showed that Mgll (A 4481) and excited 
Sill was present in the optical spectrum of this cool star. 
Since radiative levitation calculations predict that much 
higher temperatures are required before Mg can be sus- 
pended in the atmosphere, observable quantities of this 
species in EG 102 were interpreted as indicative of ongo- 
ing accretion, either from a low mass companion (for which 
no evidence was found) or from a diffuse interstellar cloud. 
Later, HBS analysed the IUE NEWSIPS spectrum of EG 
102 and noted the presence of AlII and Al III lines at the 
photospheric velocity. This observation was also interpreted 
as being a product of ongoing accretion. Photospheric and 
ISM velocities determined for EG 102 in this work are in 
good agreement with tho se measured b y HBS . 

Data provided by IChaver et alJ il995T) indicate an 
abundance of N(A1) = 10 -9 N(H) for a star with log g 
and T e fr similar to those of EG 102. The Al abundance 
in EG 102 has been estimated by fitting data with a pure 
H model spectrum into which Al has been added, using 
the code SYNSPEC. After smoothing the output spectrum 
to the resolution of IUE, abundances of order N(A1) = 
1 — > 5 x 10 -8 N(H) are obtained. The se figures are sig- 
nificantly h igher than tho se suggested bv lChaver et al.l but 
the work of IChaver et al] does not adopt the self-consistent 
appro a ch found in more recent work (e.g. iDreizler fc Schuhl 
j200lD : ISchuh fc Dreizleil (l200lT) . discussed in section 1), and 
a failure to accurately predict the Al abundance in EG 102 
cannot be used to infer the presence of accretion processes 
in the star. 

We note that lZuckerman fc Reidl il998l) find significant 
quantities of Ca in EG 102 2 (Ca/H = 2.5 x 10~ 7 , EW 
= 29mA for the A 3933 Ca K line) using high resolution 
optical echelle spectra obtained with the Keck telescope; 
current theories of radiative levitation do not predict the 
Ca ion to be present in such a cool DA. IZuckerman fc Reidl 
also find significant quantities of Ca in both of the close 
white dwarf/red dwarf pairs in their survey, and suggest 
that the presence of heavy ions such as Ca in cool DA stars 
may be attributed to binarity. To date, no companion to 



2 Listed under the alternative name G238-44 in that work. 
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Al III features in Wolf 1346 
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Figure 16. A 1854, 1862 A Al III lines in Wolf 1346, with model 
spectrum superimposed (darker line). The Al abundance inferred 
from the 1862 A line is approximately 2.2 xlO -9 N(H). Struc- 
ture in the continuum level around the A 1854 line complicated 
abundance measurements using this feature. 



EG 102 has been observed. 
Wolf 1346 

The pres ence of Si in the photosphere of Wolf 1346 was re- 
ve aled bvlBruhw eiler fc Kondol il982Tl . but later questioned 
bv IVennes et alJ <199llT who noted discrepancies between 
the velocity of Si II lines and the ground-based photospheric 
velocity, and found that abundances determined from the 
Sill lines were inconsistent with the non-detection of Si III 
in the IUE spectra of this star - observations more sug- 
ge stive of a circumstella r origin. The problem was resolved 
bv iHolberg et all Jl996h . who revised the photospheric ve- 
locity and used the advanced non-LTE code TLUSTY to de- 
rive a Si abundance 0.5 dex lower than that of lVennes et al.L 
confirming the photospheric nature of the Si lines. 

The velocities of features normally associated with the 
ISM appear to fall into two groups, at -16.19 ± 0.10 kms -1 
and -7.67 ± 3.07 kms -1 . When considering the velocities 
of isolated features, these groups appear to contain sepa- 
rate species (Ol, CII and SII at -16 kms -1 , and Sill at 
-7 kms -1 ). However, this segregation breaks down when 
those interstellar Sill features which are blended with pho- 
tospheric features, are included. For example, the ISM com- 
ponent of the A 1260.4221 Sill line, lies at -16.73 kms -1 . 

The photospheric lines (which are limited to Sill and 
Si III) show a considerable spread in velocity (from 19 to 
33 kms -1 ), with the weighted mean being v phot = 24.32 ± 
1.41 kms -1 . Given this range of velocities, the reality of 
the two ISM groupings is questionable. Data from HBS also 
show this grouping, though only four ISM lines are recorded 
(compared to ten in the current study). HBS treat these ve- 
locities as belonging to the same group, and give a weighted 
mean ISM velocity of -14.85 ± 1.50 kms -1 . For the current 
work, the analogous value is vi SM = -15.38 ± 0.95 kms - . 
Thus, both vism and v pilat are found to agree with the values 
presented by HBS. 

Lines of AlIII are clearly visible at the photospheric 
velocity. Thes e features are not n oted in previous studies 
of Wolf 1346 iHolberg et alJ dl996>) . HBS). To estimate the 
Al abundance, the A 1854, 1862 lines were reproduced by 



Star 


Species 




REJ 1738+665 


CIV, 


Si IV, 


NV, OV. 


Ton 021 


CIV, 


Si IV. 




REJ 0558-373 


CIV. 






WD 2218+706 


CIV, 


Si IV. 




Feige 24 


CIV. 






G191-B2B 


CIV, 


Si IV* 




REJ 0457-281 


CIV, 


Si IV, 


NV. 


REJ 2156-546 


CIV. 






REJ 1614-085 


CIV, 


Si IV. 




GD 659 


CIV. 






* Detected bvlHolb 


ere: et a] 


J 1200*. 





Table 3. Ionic species contributing to observed or suspected cir- 
cumstellar features. Ions responsible for principal features (i.e. 
those which are clearly identifiable without the need for co- 
addition) are indicated in bold; other species should be regarded 
as tentative identifications. 



adding quantities of Al to a spectrum generated from a pure 
H+He non-LTE model (N(He) = 10 -8 N(H)), and the out- 
put smoothed to the resolution of IUE (~ 0.2 A FWHM). 
An abundance of N(A1) = 2.2 x 10 -9 N(H) was found to 
match the observed 1862 A line, and is clo se to the value of 
~lx 10 -9 N(H) implied by the results of IChaver et all A 
greater abundance (N(A1) = 6.0 x 10 -9 N(H)) was required 
to match the 1854 A feature; however, the region around this 
line exhibits unusual structure, possibly due to instrumental 
effects. The models used to determine these abundances are 
not stratified, and hence these values are likely to be revised 
when suitable stratified models become available. The Al 
lines and synthetic spectrum are illustrated in figure [161 



3 DISCUSSION 

Of the twenty three stars considered in this survey, four were 
previously known to possess features from highly ionised 
species at non-photospheric velocities: Feige 24, REJ 0457- 
281, G191-B2B and REJ 1614-085. Four new DA white 
dwarfs may now be added to the list of those exhibiting un- 
ambiguous, highly ionised components at non-photospheric 
velocities: REJ 1738+665, Ton 021, REJ 0558-373, and WD 
2218+706. A fifth object, REJ 2156-546, shows features 
which may also be interpreted as non-photospheric, although 
data of improved S/N are required to confirm this result. A 
weak blueshifted component in GD 659 is also suggested, 
though this is exceedingly faint, comparable with the struc- 
ture of adjacent noise, and hence regarded as a tenuous iden- 
tification until further data become available. Table |3] lists 
those ions appearing as circumstellar features. Observations 
of these features are restricted to resonance transitions, and 
are most common in the CIV AA1548.202, 1550.774 doublet 
(the resonance lines are most strongly coupled to the stellar 
radiation field, and hence are more susceptible to radiative 
levitation). The velocities of interstellar, photospheric and 
non-photospheric features (if detected) have been presented 
in tabled 

Of the eleven IUE spectra considered in this study, 
only one (REJ 0457-281) shows signs of highly ionised, non- 
photospheric components. In contrast, six out of the twelve 
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Location of sample stars over nH distribution (galactic co-ordinates) 




Figure 17. The distribution of survey stars in galactic coordi- 
nates. Circles indicate stars in which circumstellar features are 
observed, with all other objects plotted as squares. Symbol size 
indicates distance ( larger = closer). Ob j ects a re plotted over the 
galactic H I map of lDickev Lockmanl ll99(t) (units are cm -2 ). 



available STIS spectra reveal such features (with a further 
detection in one of the three available GHRS spectra). The 
low number of detections in IUE data is unsurprising given 
the low signal-to-noise ratio of the instrument, and the fact 
that its resolution (0.08 A at 1400 A) is equivalent to a ve- 
locity of 17 kms -1 (compared to ~ 3.2 kms -1 with the STIS 
E140M grating, and ~ 1.28 kms -1 in the E140H configura- 
tion). The non-detection of circumstellar components in IUE 
data therefore places an upper limit on the velocity of any 
shifted features, rather than proving absence. STIS data 
may yet reveal non-photospheric features in these objects. 

The need for consistently high resolution data with ad- 
equate S/N, covering all stars in follow-up studies is clear. 
Only the highest resolution STIS data were able to show 
secondary features in the ISM lines of G191-B2B, revealing 
a possible connection between the circumstellar components 
and the local ISM. These features were not resolved in the 
STIS E140M data. However, all but one of the STIS spec- 
tra included in this survey were acquired with the E140M 
grating (PG1123+189 was observed in the higher resolution 
E140H mode, but the data do not cover the C IV or Si IV res- 
onance lines). Hence, more stars in the current sample may 
possess highly ionised non-photospheric features, at velocity 
differentials too small for the existing data to resolve. Al- 
ternatively, some stars may possess very weak circumstellar 
features which are hidden within the noise of existing data. 
Only when STIS data (ideally from the E140H configuration 
and covering the important resonance lines) are available for 
all stars in this sample, will a more comprehensive study of 
the distribution of circumstellar features in white dwarfs be 
possible. 



3.1 Influence of position and intervening ISM 
column 

3.1.1 Spatial distribution 

The distribution of survey stars is depicted in figure IT7I no 
correlation between the presence of circumstellar features 
and the position of objects on the sky is apparent. Since the 
sample stars encompass a relatively wide range of distances 
(between approximately 14 and 436 pc), the absence of such 
a positional dependence is not surprising. 



3.1.2 ISM N (HI) column density 

Line-of-sight ISM N(Hl) column densities may be estimated 
by fitting models to the observed ISM Ly-a profile (after re- 
moving the stellar contribution to l ine+continuum using a 
stell ar model), as d emonst rated by iBarstow et al] (Il994al) 
and IHolberg et al] (Il999al) . Alternatively, the H I column 
density may be determined from EUVE data using the 
strong continuum absorption below the 912 A Lyman edge. 
In this case, the column density is included as a free parame- 
ter in model spectra, whi ch are fitted to the data using y 2 re- 
ductio n techniques (e.g. IBarstow et aljll99St IHolberg et al] 
Il999al) . These methods provide the average column density, 
and are insensitive to "clumping" along the line of sight. 
As an extension to this simplification, it may be assumed 
that, ceteris paribus, a more distant star along the same or 
similar line of sight will be observed through an intervening 
column, N(HI), of greater density. The figure of interest is 
therefore the volume density n(Hl) along the ISM column, 
where, in the general case, 

N (HI) = f n(Rl)ds , (1) 
Jo 

and where s is the distance to the star. For a homogeneous 
column, n(Hl) « 7V(Hl)/s. 

Previously determined column densities are available 
for many of the sample stars, while approxi mations may be 
made for others using the synthesis maps of lFrisc h fc York! 
]l983l) . and the contour maps of IParescd l^)84l) ~ These 
columns are listed in table |1J showing that a relationship 
between circumstellar features, and the average density of 
interstellar material along the line of sight, is not observed. 
For example, both REJ 1738+665 and REJ 0457-281 show 
circumstellar features, yet their line of sight ISM volume 
densities are significantly lower than objects without such 
features. 

This result is not surprising. iDupree fc Raymond! 

<1983h find that for a DA white dwarf with T efr 
= 60,000K and log g = 8.0, the Stromgren radius 
ranges from 0.07 pc (for n(H) = 10 2 cm~ 3 ) to 30.8 pc 
(n(H)=0.01 cm" 3 ). Although their work is now dated (the 
heavy element features in white dwarf spectra are attributed 
to the ionisation of circumstellar material), these figures still 
provide a useful order of magnitude estimate for the sphere 
of influence of the white dwarf. It could also be argued that 
the quoted Stromgren radii represent upper limits, since the 
extra opacity from photospheric heavy elements would be 
expected to reduce the intensity of the radiation field at spe- 
cific wavelengths. In either case, the Stromgren radius is typ- 
ically a small fraction of the distance to the star, and since 
the distribution of material along the line of sight is unlikely 
to be homogeneous, the observed average value of n(H I) may 
not reflect conditions within the Stromgren sphere. 

3.1.3 The velocity of circumstellar features and the ISM 

For a white dwarf which is ionising nearby interstellar mate- 
rial, highly ionised non-photospheric components may be ob- 
served at velocities similar to those of the intervening ISM. 
Depending on the relative velocities of star and ISM, this 
mechanism will produce both red- and blue-shifted features. 
The majority of stars possessing circumstellar compo- 
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Table 4. H I column and volume densities for the survey stars, obtained from a variety of previous studies. Data from iFrisch fc Yorkl 
il983Ti and IParesc d ll984f) should be interpreted as broad estimations only. No data available for REJ 0948+534 or EG 102. 



Star 



N(HI) (cm" 2 ) n(HI) (cm" 3 ) 



Reference 



REJ 1738+665"-* 
Ton 021* 
REJ 0558-373* 
REJ 2214-492 
REJ 0623-371 
WD 2218+706* 
Feige 24* 
REJ 2334-471 
G191-B2B* 
GD 246 
REJ 0457-281* 
PG 1123+189 
HZ 43 

REJ 1032+532 
REJ 2156-546* 
PG 1057+719 
REJ 1614-085"'* 
GD 394 
GD 153 
GD 659* 
Wolf 1346 



< 5.0 x 



1^5x 
l-t5x 
5.8 x 

5.0 x 

> 

3.25 x 
8.5 x 
2.1 + 0.1 x 
> 1.51 x 
(1.3 + 0.7) x 
1.10 x 
9.3 + 0.1 x 
4.2 x 

4.1 x 
2.75 x 
ss 5 x 
4.4 x 

6.03 x 

3.2 x 
1.10 x 



10 18 
10 20 

10 20 

10 18 
10 18 
10 21 



10 1S 
10 18 
10 18 
10 18 
10 17 
10 19 
10 17 
10 18 
10 18 
10 19 
10 20 
10 18 
10 17 
10 18 
10 18 



1.5 

1.1 



> 



< 6.7 x 10~ 3 
-» 7.5 x 10- 1 
-► 5.5 x 10" 1 
2.73 x 10~ 2 

1.67 x 10~ 2 
> 7.4 x 10" 1 

1.35 x 10" 2 
2.65 x 10" 2 
1.4 x 10~ 2 

6.8 x 10~ 3 

3.9 x 10" 4 

2.4 x 10" 2 
4.2 x 10~ 3 
1.1 x 10" 2 
1.0 x 10" 2 

2.17 x 10~ 2 
fa 1.89 

2.5 x 10~ 2 

2.68 x 10" 3 
2.0 x 10" 2 

2.55 x 10~ 2 



Barstow et al. Q994b) 
Frisch fc York ri983'> 
Frisch fc York (19831 



Wolff et al. (19981 
Wolff et al. (19981 
Frisch fc York (19831 
Wolff et al. (19981 
Wolff et al. (19981 
Barstow et al. (19991 
F^iscione^A^jtjiL (1994 1 
Barstowe^aL 
jlolber^^^^L 
Barstqwe^aj. 
Holbcr^^^aJ^ 
Holber^^^^L 
Wolff et al. (19981 
Paresce (19841 
Barstow et al. (19961 
F^uscior^^j^et^aL 1+9941 
Holberg et al. (1999b1 
Fruscione. A. et al. (19941 




'Circumstellar features observed or suggested. 
lr Total column (Ng) quoted 



nents show little agreement between v cilc and i>ism- This 
result is unsurprising, since observed ISM absorption fea- 
tures may be blends of several unresolved ISM components 
at the resolution of the current data, and hence the dom- 
inant component (typically the LISM) may not be that in 
which the star is immersed. Nevertheless, it is useful to con- 
sider the residual value of |uis M - fcirc| (figure [T51l . These 
data show that in the majority of objects possessing cir- 
cumstellar features, the difference in velocity between these 
and the primary interstellar cloud is typically 10 kms -1 or 
less (exceptions being REJ 1614-085 and GD 659). The most 
interesting cases are found in REJ 1738+665 and G191-B2B, 
where excellent agreement is found between vi S m and D circ , 
suggesting that the shifted features arise in a Stromgren 
sphere of material belonging to the primary identified ISM 
component. 



70 
60 - 
50 - 
| 40 
> 30 
20 - 
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Velocity comparisons for stars with circumstellar components 



■ ISM - circumstellar 

□ Photo spheric - circumstellar 

□ Photospheric - ISM 



mam 




More detailed investigations are required before this hy- 
pothesis can be confirmed, and a detailed analysis of the ion 
populations and spectral characteristics of such a system 
is essential. The case of G191-B2B also acts as a caution 
against over-interpretation of these results; no correlation 
between Vism and v circ would have been recognised if not fo r 
the higher resolution data discussed bv ISahu et al] lll999l) . 
Similar correlations with hitherto undetected ISM compo- 
nents may be found in future studies based on E140H grating 
spectra. Despite the relatively narrow wavelength coverage 
available with this instrument, the current results justify a 
comprehensive program of E140H white dwarf observations, 
possibly tuned to a bandpass covering the CIV resonance 
doublet, which is most frequently accompanied by shifted 
features. 



Figure 18. Comparison between photospheric, ISM and circum- 
stellar velocities in stars showing highly ionised non-photospheric 
features (Feige 24 omitted since v ckrc and v phot vary with orbital 
position, while Vjsm is stationary). 



3.2 Metallicity and mass loss 

iMacDonaldl ]l992l) discusses the interaction between the 
flow of ISM material around a white dwarf star, and the 
weak stellar wind. In this work, the rate of mass loss, M, 
fro m the white d warf is estimated using theory developed 
bv lAbbottl ]l982h . viz. 



M w 2 x 10" 



L 



0.02 



-Mi 



yr 



(2) 
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where Z is the metallicity, relative to solar abundances, of 
a star with luminosity L and mass M. 

However, the work of lAbbottl is concerned with the en- 
velopes of O- to G-type stars, and results are found to be 
most successful for OB stars. Conversely, the theory does 
not explain mass loss rates in Wolf-Rayet stars, which are 
so mewhat different in structure. Further, since the theory 
of lAbbot J is concerned with main sequence objects, the 
metallicity parameter is formulated in terms of solar abun- 
dances, and cannot be applied directly to the broad range 
of heavy element compositions exhibited in white dwarfs. 
Hence the use of equation [5] in the current work is not en- 
tirely appropriate. 

Nevertheless, it is interesting to compare the rela- 
tive mass loss rates of sample objects calculated using 
eq uati on [5] Indiv idual abundances have been calculated 
bv lBarstow et all l)200lfl . These values were determined by 
matching observational data to a synthetic spectrum calcu- 
lated using SYNSPEC, based on a model of appropriate T e g 
and log g generated by the non-LTE code TLUSTY. This 
information is available for all objects in the sample except 
PG 0948+534 (table 0. For each star, the metallicity pa- 
rameter, Z, is calculated using the expression 



Mass loss rates in sample stars 



(3) 



where ^4*,© is the abundance of the element of atomic num- 
ber z relative to hydrogen in the star and in the Sun 3 re- 
spectively. Only "metals" (elements heavier than He) are in- 
cluded. Note that although equation |31 should be evaluated 
for all elements heavier than He, only clearly identifiable 
species as presented in tabled] have been considered here. 

Stellar luminosities quoted in tabled have been calcu- 
lated from the familiar expression 



L% = W 0A{M Q -M,} L 



(4) 



where Mq represents the absolute bolometric magnitude 
of the Su n (+4.7), or the star (calculated using the stellar 
models of lWoodlll99l . 

No clear correlation is found between M and the pres- 
ence of circumstellar features, although loss rates extend to 
lower values for objects without these features (figure H9*t . 
The lowest values are found in the coolest stars, as expected 
given the dependence of equation|2]on luminosity and metal- 
licity. However, GD 246 (T efr = 53.700K) is significantly 
hotter than REJ 1614-085 (T cfr = 38,500K) despite having 
an appreciably lower calculated loss rate. While the subset 
of stars which exhibit non-photospheric components lacks 
any object with M < 3 x 10~ 15 Mq yr _1 , (compared to a 
minimum value of 1.9 x 10~ 21 Mq yr" 1 for those with no cir- 
cumstellar components), the number of objects in this sur- 
vey is insufficient to determine the authenticity of a lower 
limit to the mass loss rates in stars showing highly ionised, 
non-photospheric components. 



3 Using data provided by iGrevesse fc Sauvall 11998 
J2 z a ( z )q ~ 1-18E-3, where z=6, 7, 8, 14, 15, 16, 26 and 
28, reflecting the elements considered in table IB1 



1.0E-11 - 
1.0E-I2 - 
I.0E-13 - 
1.0E-14 - 
1.0E-15 - 
"3 1.0E-I6 - 
1 1.0E-17 - 

I.0E-19 - 
1.0E-20 - 
1.0E-21 - 



□ □ p 



Figure 19. Calculated mass loss rates for those sample stars 
containing metals, as defined in table|F] Filled squares: stars with 
no observed circumstellar features. Open squares: stars exhibiting 
circumstellar features. Objects for which no metals are present or 
for which data are unavailable have been omitted. 



3.3 Gravitational Redshift 

The apparent velocity of absorption features formed in the 
white dwarf atmosphere will be affected by the radial veloc- 
ity of the star, and by gravitational redshifting. The velocity 
change due to gravitation will be lower in features which are 
formed in material further from the stellar surface (e.g. a 
circumstellar cloud), and will be effectively zero for a cloud 
with a sufficiently large inner radius. The gravitational red- 
shift at the stellar surface therefore defines a range of ve- 
locities, with respect to the apparent photospheric value, at 
which highly ionised non-photospheric features may be at- 
tributed to material residing within the gravitational well of 
the star. 

Gravitational redshifts can be measured directly in bi- 
nary systems such as Feige 24 (iDupree fc Raymond! Il982l : 
IVennes fc Thorstensenlll 994'), but for isolated systems, the 
velocity component due to gravitation redshift, « glav , can be 
estimated using the standard formula, 



Xob 



A y 



Art 



GM 



(5) 



where G is the universal constant of gravitation, M is the 
stellar mass, R is the radius (of the star, or of the circum- 
stellar cloud), c is the speed of light, and \ bs,rest are the 
observed and rest wavelengths of the absorption line, re- 
spectively. Substituting values for the physical constants, 
the expression for w grav becomes simply 



6.36 x 10' 



M, 

"r7 



(6) 



where solar units are to be used for the white dwarf mass 
and radius. In the work which follows, estimates using this 
expression are used to identify features which could arise 
from material residing well within the gravitational potential 
of an object. 

Values for « grav , calculated using equation [5] are in- 
cluded in table H The presence of material in the gravi- 
tational potential well of a star provides no explanation for 
objects in which material appears at a velocity redshifted 
with respect to the photospheric value, as in the case of WD 
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Table 5. Photospheric heavy element abundances of the surveyed stars, determined from far-UV spectroscopy (no data available for 
PG 0948+534). Note that gaps in the table for those stars up to (and including) REJ 0457-281 are mostly due to the absence of data 
for a particular spectral range rather than a true absence of the element itself. F or the remaining, cool er stars, gaps in the table reflect 
genuine absence of these species at abundances detectable by IUE or STIS. From lBarstow et aljl200lh . 
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*Circumstellar features observed 



2218+706, Feige 24, REJ 2156-546, and REJ 1614-085. The 
mechanism may also be ruled out for three further stars, 
REJ 1738+665, Ton 021 and REJ 0457-281, in which the 
calculated gravitational redshift is substantially lower than 
required to explain the velocity of non-photospheric compo- 
nents. 

For the three remaining stars, the gravitational red- 
shift is comparable to the velocity difference between photo- 
spheric and circumstellar components. However, if the cal- 
culated values of u grav are assumed to be reasonable, only 
circumstellar features in REJ 0558-373 may be explained by 
the presence of matter inside the gravitational well (at ap- 
proximately 5 stellar radii from the surface). In the cases 
of G191-B2B and GD 659, the non-photospheric compo- 
nents differ from v v h<,t by an amount equal to, or slightly 
greater than, the value of v grav . This suggests that the ma- 
terial lies at a radius greater than that at which gravitational 
redshifting produces observable changes in line velocity. It 
is therefore apparent that, with the possible exception of 
REJ 0558-373, blueshifted features in the spectra of objects 
in this study cannot be explained by the presence of material 
within the gravitational potential well of the star. 

3.4 Non-photospheric material and its relation to 
planetary nebulae 

The existence of an old planetary nebula (PN) around two 
of the survey stars has already been addressed in the cases of 
REJ 1738+665 and WD 2218+706. Mass loss and the pro- 
duction of a PN as a white dwarf progenitor leaves the AGB, 
are relatively well accep ted processes (thou gh the precise de- 
tails are certainly not (lLangill et alJll994l) ). It is therefore 



pertinent to ask whether the presence of non-photospheric, 
highly ionised features around other DA stars at the hotter 
end of the sample (i.e. those for which the nebula material 
may not be completely dispersed), is consistent with mate- 
rial from a (now ancient) PN around the star. 

Expansion velocit ies for PN e are available in 
the literature (e.g . IWjanberger] [l9 89.1. The study 
by iNapiwotzki fc Sclioiip'eni^illu995l) is particularly 
relevant, since it deals specifically with old planetaries, 
with central stars in the advanced stages of transition 
into the white dwarf area of the H-R diagram. Since even 
the hottest stars in the curr ent study are more highly 
evolved than those covered by INapiwotzki fc Schonbernerl 
average temperatures are cooler, and any nebular material 
presumably of lower density and more widely dispersed. 
Direct detection of nebular emission around these stars 
is therefore difficult, excepting the two cases previously 
highlighted. Further, in the case of the older, cooler stars in 
the current work (e.g. REJ 2156-546 and REJ 1614-085), 
the planetary nebula must have dispersed long ago, and is 
therefore unlikely to offer an explanation for the existence 
of highly ionised, non-photospheric features. Nevertheless, 
the possibility that some circumstellar features are of this 
origin may be ass essed by comparison with th e expansion 
velocities noted bv INapiwotzki fc Schonbernerl 

Figure 12771 shows the distribution of expansion velocity 
(in this case plotted against nebula radius ) for the objects in- 
vestigated bv lNapiwotzki fc Schonbernerl No data are avail- 
able for the radius of any PNe which may surround stars 
in the current survey; however, an analogous velocity can 
be derived in the form of the difference between the veloc- 
ity of circumstellar and photospheric features, as listed in 
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Planetary nebula expansion velocity versus radius 



WD22I8+706 
(Napiwotzki & 
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Figure 20. Circles: Nebula expansion velocity v exp 
plo tted against nebula r adius, for stars studied 
by iNapiwotzki fc Schonberncr il995f) . Triangles: tt phot - v cilc 
for stars in the current work (no values for radius are available). 



absorbing matter, and the old, dispersed PN material sur- 
rounding these stars. 

The results of this comparison suggest that further work 
on the link between highly ionised non-photospheric fea- 
tures, and planetary nebula is justified. Although this hy- 
pothesis appears to contradict the suggestion that shifted 
features may be related to interstellar material within the 
Stromgren sphere of a white dwarf, both mechanisms may 
operate in different objects, with Stromgren spheres as the 
dominant source for cooler objects. Correlating typical PNe 
densities with the column densities derived from curve of 
growth analysis would provide further evidence of a relation- 
ship between these two apparently separate entities. This 
comparison is complicated by the considerable difficulty in 
estimating the masses of PNe, arising from uncertainties in 
the distances to nebulae, and from the ongoing debate as 
to whether these objects are typically ionisation or mass 
bounded. 



table Note that WD 2218+706 (DeHt5) is common to 
both studie s, and in this case the true v alue of v exp , as de- 
termined bv lNapiwotzki fc Schonberneri is considerably dif- 
ferent to the negative velocity, relative to the photospheric 
value, found in the current work. It is therefore clear that 
the non-photospheric features of WD 2218+706 discussed in 
the current study are not produced by the observed plane- 
tary nebula (though the nebula cannot be discounted as a 
source of this material). 

We note that the observation of features which 
are redshifted with respect to the photosphere is not 
inc ompatible with the nebula hypothesis; for exam- 
ple, |^weg(^^^^^iwotzkj| (Il994l) discuss the planetary neb- 
ula Sh 2-174 and the white dwarf GD 561 observed on one 
edge of this nebula. The objects are found at similar dis- 
tances, and [OIII] emission in the nebula is located imme- 
diately adjacent to the white dwarf. These observations, the 
statistical improbability of GD 561 being an isolated hot 
white dwarf which happens to be wandering through the 
nebula, and the difficulty in explaining the existence of a 
small nebula other than being of PN origin, confirm that 
GD 561 is indeed the source of the planetary nebula. The 
distinctly non-spherical morpho logy of Sh 2-174 is far from 
unique in studies of old PNe (iTweedv fc Napiwotzkil and 
references therein). Gross asymmetries are attributed to in- 
teraction between the nebular material and the surrounding 
ISM through which it moves. Given the location of GD 561 
relative to the central region of Sh 2-174, it is clear that 
the PN material may produce spectral features which are 
either blue- or red-shifted with respect to the photosphere, 
depending on the angle from which the system is viewed. 

For all but two of the stars in this work, v vtlot - v cilc is of 
a similar order of magnitude as the expansion velocity typ- 
ical of old PNe. Two stars (WD 2218+706 and REJ 2156- 
546) are shown with the negative values of v ciTC discussed 
above - but the absolute values, |« P h 0t - v cilc \, are consistent 
with the remainder of the sample. In contrast, radiatively 
driven winds from the surface of a white dwarf should be of 
a similar ord er to the stellar es cape velocity, i.e. ~ 1000 km 
s or more iMacDonaIdlll992l) . These results may indicate 
some form of link between the origin of the non-photospheric 



4 SUMMARY 

We have described the detection and interpretation of highly 
ionised absorption features at non-photospheric velocities, in 
high resolution UV spectra of hot DA white dwarfs. These 
features may be indicative of accretion or mass loss in white 
dwarfs - processes which may explain the non-equilibrium 
abundances, compared with the predictions of radiative lev- 
itation theory, observed in many objects. Four of the stars in 
the sample were previously known to show non-photospheric 
features: Feige 24, REJ 0457-281, G191-B2B and REJ 1614- 
085. This work has revealed at least four new objects with 
multiple components in one or more of the principal reso- 
nance lines: REJ 1738+665, Ton 021, REJ 0558-373 and WD 
2218+706. A fifth object, REJ 2156-546 also shows some ev- 
idence of multiple components, though further observations 
will be required for their reality to be confirmed. 

Several possible mechanisms for the formation of these 
features have been discussed. The presence of material 
within the gravitational potential well of a white dwarf is 
found to be an unsatisfactory explanation for the produc- 
tion of these features. Predicted mass loss rates based on 
the luminosity and metallicity of stars show no correlation 
with the presence of shifted features. However, these mass 
loss rates are calculated using theories developed for main 
sequence stars, and may be inappropriate for application to 
highly evolved objects. Further, the quantification of metal- 
licity is a highly subjective measurement, and is likely to be 
a major source of uncertainty in these calculations. 

A possible correlation is observed between the velocity 
of shifted features and that of the ISM. This is particularly 
obvious in REJ 1738+665 and G191-B2B, which show very 
close matches between Vism and u circ . For most of the re- 
maining stars, the difference between these velocities is less 
than 10 kms -1 . Higher resolution observations are required 
to detect the presence of multiple ISM velocity components, 
which may reveal further correlations - as demonstrated by 
the case of G191-B2B. 

An alternative or additional source of shifted features 
may be found in planetary nebulae. Velocities of shifted fea- 
tures with respect to the photosphere in this study are found 
to be entirely consistent with the expansion velocities typ- 
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ical of old PNe. By appealing to the irregular morphology 
of highly evolved nebula, both blueshifted and redshifted 
features may be explained. Detailed modeling of the inter- 
action between the white dwarf and surrounding material 
should determine whether stellar radiation alone is suffi- 
cient to produce the observed ionisation, or whether addi- 
tional excitation (perhaps in the form of shock-heating) is 
required (Napiwotzki, private communication). 

The non-detection of highly ionised non-photospheric 
features in many of the stars investigated may indicate their 
absence, but equally, may reflect the limited resolution and 
signal-to-noise ratio of available data. This is particularly 
important when considering non-detections in IUE data, 
where velocity differentials of less than 17 kms -1 between 
photospheric and shifted components are below the resolu- 
tion limit of the instrument, and where the S/N ratio is 
inferior to that of more modern instruments such as STIS. 
This highlights the importance of acquiring consistently high 
resolution data for all stars in this and future samples. Four 
(and possibly five) new identifications of circumstellar fea- 
tures have been made using medium resolution (E140M) 
STIS spectra; these successes demonstrate the value in high 
resolution studies of this type, and justify an extension of 
the program to include higher resolution data of improved 
S/N for an expanded sample of objects. 
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